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We consider the pricing of European-style structured credit payoff 
in a static framework, where the underlying default times are indepen- 
dent given a coinmon factor. A practical application would consist of 
the pricing of nth-to-default baskets under the Gaussian copula model 
(GCM). We provide necessary and sufficient conditions so that the 
corresponding asset prices are martingales and introduce the concept 
of "break-even" correlation matrix. When no sudden jump-to-default 
events occur, we show that the perfect replication of these payoffs un- 
der the GCM is obtained if and only if the underlying single name 
credit spreads follow a particular family of dynamics. We calculate 
the corresponding break-even correlations and we exhibit a class of 
Merton-style models that are consistent with this result. We explain 
why the GCM does not have a lot of competitors among the class of 
one-period static models, except perhaps the Clayton copula. 
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1 Introduction 



The risk management of structured credit products has become a key issue 
for many financial institutions. Even if their payoffs are only driven by the 
realization of default events, the value of these products are exposed to credit 
spread volatility. The recent credit crisis has highlighted the need for pricing 
models which can produce dynamic hedging strategies that would replicate 
the payoff at maturity of these products and, more importantly, reduce to 
zero the "P&L volatility" of structured credit trading portfolios. 

The most common structured credit products are synthetic CDO tranches 
(also known as CSO tranches). A standard modelling approach to value CSO 
tranches is the one factor Gaussian copula model (GCM). This model spec- 
ifies directly the joint law of the underlying default times (Li, 2000), where 
the dependence between these default times is measured through a corre- 
lation parameter. With the marking of correlation levels of base tranches 
[0,i^'], K G [0,100%], GCM has become a market standard known as the 
"base correlation" model (see O'Kane and Livasey 2004, for instance). 

The Gaussian copula approach has been criticized for several reasons. 
Firstly, the model does not define any credit spread or default intensity dy- 
namics. It only takes as a set of inputs current credit spread curves, recover- 
ies and a correlation. Secondly, the common opinion is that CDO tranches 
are not priced through a replication argument, contrary to the Black-Scholes 
model. D. Li also wrote: "The current copula framework gains its popular- 
ity owing to its simplicity. However, there is little theoretical justification of 
the current framework from financial economics. We essentially have a credit 
portfolio model without solid credit portfolio theory". In our opinion, these 
criticisms are partly due to a lack of understanding of the properties and 
limitations of the GCM. In this article, we will consider the larger class of 
static factor models, for which the default times are assumed as independent 
conditionally on a random variable (the so-called "common factor"). We will 
explain to what extent and under which conditions such copula models can 
be seen has pricing models by replication, as in the standard Black-Scholes 
theory. In a continuous spread model, we show that the payoff of an arbitrary 
European basket credit derivative can be perfectly hedged under particular 
spread dynamics and if the copula pricing model satisfies some functional 
equation. 

The replication of CDO-type cash flows has generated a significant amount 
^in The Definitive Guide to CDOs, G. Meisner (ed), Risk Books, p. 71 
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of literature in the last years. Some authors have tried to write explicit hedg- 
ing strategies for some particular top-down models, particularly when risks 
are driven mainly by sudden default events (pure jump dynamics with con- 
tagion): Arnsdorf and Halperin (2008), Frey and Backhaus (2008, 2009), 
Herbertsson (2008), Laurent et al. (2008), among others. Then, the relevant 
hedging instruments are most often instantaneous CDS (virtual instruments) 
and/or credit indices. See the survey of Cousin, Jeanblanc and Laurent 
(2010). Recently, Cousin and Jeanblanc (2010) have obtained general the- 
oretical results when the default events are ordered. None of these authors 
deals with the case of static factor models, particularly the market standard 
GCM. Others have studied empirically the performances of various delta 
hedging strategies, including the current practice: Cont and Kan (2008), 
Ammann and Brommundt (2009), Cont et al. (2010), Cousin, Crepey and 
Kan (2010) etc. In this article, we choose rather an opposite point of view: 
we restrict ourselves on the diffusive part of the risk (spread risk), but under 
standard market models and with individual CDS as hedging instruments. 

In section|2] we specify the framework and the notations. We explain how 
the price given by a general factor model can become a martingale. In sec- 
tion [3j we particularize the latter result in the case of the one- factor Gaussian 
copula model (IFGCM). We exhibit the single credit spread dynamics that 
are consistent with the pricing by replication of credit basket structures in 
this model. We call "break-even" correlation the flat Gaussian copula corre- 
lation level to perfectly replicate a given structured product with individual 
CDS, when no sudden default occurs. Break-even correlations are written 
explicitly as some functions of the spread volatilities and correlations. Sec- 
tion [1] extends these results in the general case of p-factor Gaussian copulas, 
and introduces the concept of break-even correlation matrix. In section [5} we 
calculate explicit hedging errors for all First-p-To-Default securities. In sec- 
tion [6j we show the one-to-one correspondence between our spreads dynamics 
and a class of structural models . Alternative factor models are studied in 
section [7] Finally, section [8] provides empirical illustrations of our results, 
particularly break-even correlation levels, by simulating spread trajectories 
and by using real market quotes. 

2 The framework 

First, let us describe our underlying market and technical assumptions. To 
simplify, we will assume we are not exposed to any interest rate risk: 

Assumption (I): interest rates are zero now and in the future. 
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The basic building blocks in the market are n individual T-maturity 
Credit Default Swaps, or, equivalently, n survival probabilities up to T. To be 
specific, for every name i and every maturity T, let the survival probabilities 
Qit{T) be 

Qu{T) = Qin > T\Tt) := exp {-hu.{T - t)) . (1) 

When there will be no ambiguity, such survival probabilities will be denoted 
by Qit or even Qj. In the equation above, Q denotes a "risk-neutral" proba- 
bility measure and hu is the hazard rate of i between time t and maturity T. 
From the quantity Qu, we can deduce the price of the T-maturity CDS of 
name i: for every i and t, dCDSu = —dQit{T).{l — Ri), where Ri is the con- 
stant recovery rate of i. As a result, Qu can be considered as a traded asset 
and our market will be defined as (Qit{T), . . . , Qni(7'))tg[o,r]) endowed with 
its natural filtration := {J''t)te[o,T]j -^t = o'iQiu{T),u < t,i = l,...,n). 
Note that the filtration T contains all the relevant information at time t, 
including past default events as well as past and current credit spreads. 

Second, let's consider a credit basket product of n underlying names, 
n > 2. Its payoff depends only on the realizations of default events in the 
underlying basket up to the maturity T of the basket. We will further assume 
that all the default payments are made at maturity and that the premium to 
enter this derivative contract is paid upfront. Mathematically, we can write 
the payoff using indicator functions: 

i^{Si{T),...,5n{T)), 5i{T) = l{Ti<T), i = l,...,n, 

where Tj is the default time of obligor i. For convenience, set 5{t) = {Si{t), . . . , 6n{t)). 
This framework encompasses potentially all the European credit basket deriva- 
tives, through an Arrow-Debreu type payoff decomposition. 

Placing ourselves within the framework of an arbitrage- free financial mar- 
ket model, as in Bielecki and Rutkowski (2002) for instance, the t-value of 
any attainable product is the expected value of its future cash-flows under a 
risk-neutral measure Q. Here, it will be 

Vt = ^mT))Q{6{T)\Tt). (2) 

S{T) 

In other words, (Vt)(g[o_T] has to be a J^-martingale under Q. 

In order to derive tractable pricing and hedging of the basket derivatives, 
we make a further assumption that will allow us to work in the class of "static" 
factor models: 
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Assumption (Fl): At any time, conditionally on some random factor 
X £ MP and the current market information, the underlying default events 
(tj < r), i = 1, . . . , n, are independent. The factor X has got a density fx 
w.r.t. the Lebesgue measure in M^. 

With the latter assumption, we encompass most of the current market 
models in the credit area (see for instance Burtschell et al. 2005). Note 
that "systemic" factor X depends on the current time t implicitly, but we 
specify no dynamics concerning the successive X. This factor X could be 
multidimensional and when X is a vector, its values will be denoted in bold. 
When we restrict ourselves to a one factor model, we will add this to our set 
of assumptions: 

Assumption (F2): The dimension of systemic factor X is one. 

Under (Fl), the risk-neutral joint law of the default times is defined 
by some intermediary quantities called "conditional default probabilities" 
Pi\x = QiTi < T\X = x,J^t), that depend on Qn, t and x only. Similarly 
the "conditional survival probabilities" will be denoted by q^^ := 1 — p^^. 
It is proved in section [g] that (Pi|x)te[o,T] is a martingale under Q and for 
a convenient extension of the filtration J^. There, x will be the proper t 
realization of an adapted process in this extended filtration. 

The t- value of the structured product is denoted by Vj. As usual in a 
factor model and under (Fl), 

„ n 

S(T) J=l 

where the summation is w.r.t. all the n-vectors of default indicators. Note 
that Vt = V{Qit, ■ . . ,Qnt), and Vt does not depend on time or remaining 
maturity, except through the survival probabilities themselves. Equivalently, 
there is no theta effect: dVt/dt = 0. 

The key issue we are interested in is the following one: Under which 
conditions is the "ad-hoc" standard pricing formula ([3| a J^-martingale ? In 

other words, can we write Vt = Eq 

claim (see (121))? 



^{6)\J-'t , as we can for any attainable 



Note that the pricing formula in equation (|3]) does not depend on the 
volatilities of the default intensities and/or probabilities, due to the static 
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nature of the underlying model It does not mean that the volatilities 
Git do not matter to calculate the price Vt- Indeed, we will see that they 
influence the choice of the "right" pricing parameters, particularly the pricing 
correlation in the GCM. 

Assumption (A): no default event will occur in the time interval [0,T). 

Apparently, the latter assumption is a bit provocative: it assumes we are 
living in a credit world without credit events except at the maturity T. How- 
ever, this does not preclude that the spread of any name becomes so wide 
that it trades close to default within that period. Also, for short/medium 
term horizons and for investment grade portfolios (like Master indices), the 
assumption (A) will be acceptable and is equivalent to consider that the 
structured product, while in a trading book, is exposed to spread risks 
mainly. It could be stressed that, historically, most of the default events 
have been anticipated by the main market dealers some weeks before their 
official announcements with the price of the CDS trading near the recovery 
value. Additionally, we assume that the survival probabilities Qit(T) are Ito 
processes. Thus, we will live in a Black-Scholes world. 

Assumption (Q): Under the historical measure P, the survival probability 
process of name i up to maturity T is given by dQit{T) = fJ-i^t,Tdt+ai^t,TdWit, 
when i is not defaulted, and Qit{T) = else. The volatility and drift pro- 
cesses above are J^-adapted. Obviously, the processes {Wn) are correlated 
J^-Brownian motions under P: Ep[dWit-dWjt] = Pijdt, i ^ j. We will call 
Pij the "spread correlation" between the names i and j. 

As deduced from standard arguments, the risk neutral Q can be assumed 
to be the risk-neutral measure given by the Girsanov's theorem. Therefore, 
under Q, the survival probabilities follow the processes dQit{T) = cri^t,TdWit, 
when i is not defaulted, and Qit{T) = else. Here, the processes {Wit) are 
correlated J^-Brownian motions under Q, and EqldWit.dWjt] = Pijdt, i 7^ j. 
Since there will be no ambiguity, (Ti^t,T is denoted by an simply. We denote 

'■= [Pij]i<i,j<n, where pu = 1. We prove 



^This is different from the Black-Scholes model, where the pricing of options involves 
some volatility parameters. Potentially, the latter ones are different from the "true" realized 
volatilities, inducing some hedging errors: see Carr and Madan (2001). In our case, such a 
discrepancy between the risk neutral parameters and the realized ones will appear below, 
but through correlation levels : see section [Sl 
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Theorem. 1 Under the assumptions (I) (Q) (A) (Fl) (F2), 



"^CTitcrjtPij 



5{T) 

dQi dQj 



{i,j),i<j 



fxix) 



(7, 



' dx 



Xiix) 



2 dpi\a 

dx 



Xj{x) 



dx + rj 



where 



and 



dt 



2^ 

5{T) 



fx{u) du, 



n 



1-5, 



Pj\x1j\x 



+00 



When there are two names only in the basket, the product over k above is 
just replaced by the constant one. See the proof in the appendix [A} The last 
term residual term rj is zero in a lot of models, particularly the Gaussian 
copula one. Actually, rj = when lim|j.|_j._|_oo Xiix) = for every name i. 

Therefore, under the "no default" assumption (A), it is possible to get a 
strictly flat P&L by delta hedging any structured product with T-maturity 
individual CDS if and only if the following partial differential equation is 
satisfied: for every couple of indices (i, j), i ^ j, for every time t and every 

X, 

Note that, in general, traders try to be hedged simultaneously against 
spread moves and against jump-to-defaults. This means they have to trade 
two types of hedging instruments, typically two CDS per name with different 
maturities. In that case, it is no more possible to get a flat P&L, because the 
trader will pay carry to be protected against default events. In average and 
in the long-term, this additional cost is zero, but not between two successive 
times. 



(4) 



3 The Gaussian Copula revisited 

It is now fruitful to consider the Gaussian copula model under the framework 
of section [2] Indeed, it is still the market standard, even if it has been 
criticized, particularly during the credit crisis 2007-2009. 
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Assumption (Cop): To price our European structured product, the frame- 
work is the one-factor Gaussian copula model. The correlation between the 
i-th asset value and the systemic factor X is denoted by pi. 

Since we consider payments at maturity only, the timing of default events 
and term structures of credit spreads do not matter to price the structured 
product. As a consequence, (Cop) is equivalent to define the joint law of 
the default events 6{T) only. In general, both approaches are not equivalent. 
Typically, the so-called "base correlation" methodology specifies the joint law 
of default events at several maturities, but it is not possible to infer from 
them the joint law of default times, except when the correlation levels are 
the same for all maturities. 

As a consequence of (Cop), the correlation between two different asset 
values Xi and Xj will be piPj. Set ^ = 1 — <I> and 

Under (Cop), the "conditional default probability" of name j by time T is 
by definition 

and qj\^ = I - Pj\x = ^{dj). 

It is proved in the appendix [B] that the function Xj of theorem [T] is now 

^^^^^-</.(ci>-i(g,,)) dQ, 

Equation ([6| is a noteworthy and a nice property of the Gaussian copula 
specification. Note that Xji^) tends to zero when \x\ tends to the infinity. 
We deduce that the "noisy" term r] in theorem [T] is zero. 

Theorem. 2 Under the assumptions (Cop) (I) (A) (Q) (Fl), 

dv, = i...).dm + ^Yl E / n <4''-(2<^^ - i)(2'J. - 1) 

S{T) «<i 



Pj\l^- PI 
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where 



Pit : 



In practice, it is highly valuable to exhibit the "right" values of pricing 
parameters, i.e. the values under which (Vt) becomes a J^-martingale under 
Q. Then, without any default default before T, we can hedge our struc- 
tured product dynamically with T-maturity CDSs'. Under (Cop), the single 
parameters are the "beta" factors pi, i = 1, . . . ,n. Remind that we know 
the "real- world" correlation parameters pij. This leads us to introduce the 
concept of break-even correlation. 

Definition. 3 A given structured credit security is priced under the one- 
factor Gaussian copula framework with a constant parameter p. This security 
is hedged dynamically with individual T-maturity CDS and no default occurs 
during the time interval we consider. Then, the value of the copula (fiat) 
correlation that induces a zero P&L in this time interval and for a realized 
trajectory of individual credit spreads is called a "break-even" correlation. 
When the time interval is infinitesimal, we will be calling it "instantaneous 
break-even" correlation. When its value is independent of the (random) credit 
spread trajectories, it will be called a "universal break-even" correlation. 

For us, the term "break-even" correlation is related to p'^^ ^ and not to 
PBE,t- The latter quantity will be called "beta factor" rather. This concept 
has been introduced first in |5j for First- To-Default. A break-even correlation 
is comparable to a base correlation. Its squared root is the so-called "beta"- 
factor used by practitioners, i.e. the correlation between an individual asset 
value and the systemic random factor in the one-factor Gaussian copula. 

Instead of looking for a single break-even correlation number and flat cor- 
relation structures, we can extend these definitions straightforwardly to deal 
with one-factor pricing structures. Then, the beta factors pi, i = 1, . . . ,n, 
that induce a zero P&L in a given time interval and for a realized trajectory 
of individual credit spreads will be called "break-even beta factors". 

Invoking theorem 2 the break-even correlation p'^^ j in the infinitesimal 
time interval [t, t + dt\has to be a linear combination of the pairwise spread 
correlations: 
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where the weights Wij^t involve individual survival probabilities at t, instan- 
taneous volatilities an and the payoff functional |^ The weights Wij^t could 
be negative because of the default indicator functions 6k- Thus, the exis- 
tence of pBE,t is not guaranteed. The only way to be sure that there exist 
a break-even correlation, or at least break-even beta factors is to cancel all 
the terms 

2l3itl3jtpij - pipjPft - Pipjf3jt, i + j, 

for every /3it and f3jt. If it can be done, we find universal break-even beta 
factors. And they will be the same whatever the payoff because the Pa 
terms depend on the spread dynamics only. This intuition will be formalized 
in the next theorem. Before stating this result, we need to define the names 
that have really an influence on the realized payoffs of the structured product. 

Assumption (P): for every name i in the basket, there exists another 
index j ^ i and some indicators 5^, k ^ i,j, such that 

V'Cl, 1, <r* (,^^.)) - ^P{1, 0, - V(0, 1, + ^(0, 0, 5* (,,^.)) / 0. (7) 

In the equation above, the first (resp. second) index is related to name i 
(resp. j), and <J_(ij) is the n — 2-vector of indicators 5^, k ^ 

The equation ([7| means: the knowledge of i's default is informative to 
evaluate the change of the final payoff induced by j's behavior, under some 
scenario of other default events. Such a technical condition is satisfied in 
practice by all the usual credit derivatives, particularly First-To-Default se- 
curities. Since every structured credit product can be seen as a linear com- 
bination of First-To-Default on some sub-portfolios (see Brasch 2006), (P) 
is satisfied generally when i and j take really part of the payoff ■0. 

Moreover, to state the necessity part in the next theorem, we need an- 
other technical condition. 

Assumption (B): for every name i in the basket, the function Qi i— )• 
Tit{Qi) / (Qi)) from (0, 1) to M"*" is bounded from above. 

The latter assumption avoids some pathological behaviors, when survival 
probabilities become close to zero or one. It provides a type of upper bound 
in terms of volatility explosion in the tails. 

Now, we can state a striking result: we are able to exhibit the dynamics 
of Qit (under the historical and risk-neutral measures), that allow perfect 
replication in our framework. In section |C] we prove the following result. 

^Actually, the break-even correlation itself is hidden in Wij through di and dj. This 
means we have to solve an implicit function to calculate break-even correlations. 
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Theorem. 4 Under the assumptions (Cop) (I) (Q) (Fl) (A) (P) (B), any 
credit basket derivative can he perfectly hedged against credit spread variations 
if and only if 

(i) the dynamics of the survival probabilities under P are given by 

dQit = aiCt(l){^'HQit))dWit + fiitdt, (8) 

for every index i. Above, fi, denotes some positive constant and ^ a 
deterministic function such that Ymit^T /q '^'^ = +00. 

(ii) there exist some real numbers pf such that 



Pij = Ep[dWit.dWjt]/dt 



2aiO'j 



(9) 



(Hi) for every index i, the correlations that are used for the pricing of the 
structured product, in the one-factor Gaussian copula model, are given 
by the relations pi = pf . 

Corollary. 5 Under the assumptions of theorem^ assume that the matrix 
of spread correlations S"^ is one-factor, i. e. there exist coefficients pf such 
that Pij = pf Pj for all i ^ j- Then, if the coefficients an are the same for 
every index i, the pricing by replication under the IFGCM is possible. In 
this case, pi 



In general, an arbitrary correlation matrix S"^ (even one-factor) and ar- 
bitrary deterministic "volatility" coefficients cjj do not imply the existence 
of "break-even" pricing coefficients pi. Indeed, we are not sure that there 
exist coefficients pf that satisfy the key equation In other words, the 
correlation matrix T,^ must have got a particular structure, given by so 
that there exist break-even beta factors. This family of correlation matrices 
is parameterized by the "volatility" coefficients cjj, except when all of them 
are equal. 

Actually, consider the matrix S"^, obtained by reverting equation i.e. 



Pij 



2aiaj 



Fortunately, T,^ is a correlation matrix. This is proved incidentally in sec- 
tion [6] Nonetheless, we are not insured that this correlation matrix is one- 
factor. Thus, the existence of the coefficients pf is not guaranteed. 
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In section [6j we will see that, typically, the classical one-period structural 
model induces the choice = {T — t)~^^'^ and = 1. In this case, the 
only discrepancy between the underlying spread dynamics is due to different 
survival probabilities Qu and different spread correlations. Now, we have 
stated that perfect replication may be compatible with an additional level 
of heterogeneity between spread volatilities, i.e. when the are different 
among the underlying names. For instance, all other things being equal, 
when a "volatility" factor ai goes up (resp. down), then the "break-even" 
beta factor of i goes down (resp. up), so that equation ^ remains satisfied. 

Nicely, equation ^ can be solved explicitly. Indeed, it can be checked 
that the solution of the latter equation is 



Q,, = $ |^exp(^ii^)ci>-i(Q.o) + j[*exp(^^i^){a,e.dT4^» 
' du}), (10) 



where (Wu) is a Brownian motion under P. Note that the consistency condi- 
tion liuit^T Qit{T) € {0, 1} is satisfied through the property limt^T Jq £,u du = 
+00. Finally, note that the previous individual dynamics can be rewritten 
easily in terms of hazard rates (or credit spreads, approximately) hu, defined 
by Qit ■■= exp{-hit.{T - t)). 



4 Extension to the j9-factor Gaussian copula model 

Now, we extend the results of the previous section to deal with more general 
Gaussian copulas, i.e. not only one-factor. By definition, this means that 
the default events are independent w.r.t. a p-dimensional standard Gaussian 
random variable X, and that there exist vectors pi S M^, i = 1, . . . ,n, such 
that "conditional default probabilities" are 

P.|x = Qin < T\X = X, J-i) = CD ^ ^"J^g^^^ ^ := $(d,), (11) 

for alH = 1, . . . , n and x G W, and with our previous notations. In the equa- 
tion above, \pi\ denotes the Euclidian norm of the p- vector \pi\. Obviously, 
joint default probabilities are obtained by an integration w.r.t. the factor 
density fx, where /x : — >• M+, x i— )• Y\^^i4>ixk)- In this framework. 
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the pricing parameters are the vectors pi G MP, i = 1, . . . ,n. They allow to 
build a correlation matrix S of rank p, for pricing purpose. It is defined by 
Sjj = p[pj, i ^ j and Tin = 1. Thus, in this section, we assume 

Assumption (pCop): To price any European structured product, the 
framework is the p-factor Gaussian copula model, for some p G {1, . . . , n}. 

When p > 2, we cannot invoke theorem [T] anymore. Indeed, it is based on 
a single integration by parts argument. When x is replaced by a vector x, it 
is far from obvious to exhibit the relevant univariate variable of integration, 
and then to find the new functions Xi- Nonetheless and nicely, we prove here 
that the results of theorems |2] and |4] are still available. 

Theorem. 6 Under the assumptions (pCop) (I) (A) (Q) (Fl), 

• [2PuP,,p., - f/,p,0i - p'^p.^l] ^(jO^^.O^W ^ ^^^^ 
where 

See the proof in appendix [D] Assume we know the "real- world" correlation 
parameters pij and the volatility parameters This leads to the concept 
of break-even correlation matrix. 

Definition. 7 A given structured credit security is priced under the p-factor 
Gaussian copula framework. This security is hedged dynamically with indi- 
vidual T-maturity CDS and no default occurs during the time interval we 
consider. Then, the p-factor correlation matrix that is used for pricing pur- 
pose and that induces a fiat P&L in this time interval and for a realized 
trajectory of individual credit spreads is called a "break-even" correlation ma- 
trix. When its value is independent of the ( random ) credit spread trajectories, 
it will be called a "universal break-even" correlation matrix. 

Note that this matrix does not always exist for arbitrary spread correla- 
tions Pij and arbitrary volatilities au- As in the univariate case, we can detail 
necessary and sufficient conditions for the existence of a break-even correla- 
tion matrix. With the same reasoning as in theorem [4j we find the single 
dynamics that are compatible with perfect replication of credit structured 
products. 
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Theorem. 8 Under the assumptions (pCop) (I) (Q) (Fl) (A) (P) (B), con- 
sider a particular credit basket derivative. It can be perfectly hedged against 
credit spread variations if and only if 

(i) the dynamics of survival probabilities under P are given by 

dQit = ai^tH^-HQit))dWu + ^l^tdt, (12) 

for every index i. Above, fij denotes a positive constant and ^ some 
deterministic function such that lim^^^j^ du = +00. 

(a) there exist some vectors in W such that 

(pfypl yi^j. (13) 

(Hi) for every index i, the correlations that are used for the pricing of the 
structured product, in the p-factor Gaussian copula model, are given by 
the relations pi = pf , i = 1, . . . ,n. 

Corollary. 9 Under the assumptions of theorem^ assume that the matrix 
of spread correlations T,^ is p-factor, i.e. there exist p- dimensional vectors 
pf such that pij = {pf YPj for all i j. Then, if the coefficients (Ti are the 
same for every index i, the pricing by replication under (pCop) is possible. 
In this case, pi = pf . 

The result above can be rewritten in a more striking way. If the ma- 
trix T,^ = [{ -2'f-2 }Pij\ is given by a p-factor correlation structure, then the 

Gaussian copula pricing model should be p-factor. And the corresponding 
pricing correlations are given by the decomposition of S'^, that is a corre- 
lation matrix (see section [6]) . This is clearly in contrast with the current 
practice, where practitioners use only one-factor Gaussian copula model by 
far. 

5 Analysis of a First ^-th to default 

The simplest basket credit derivative is a First- To-Default, where the payoff 
is non zero, except when all the names have survived until the maturity date 



Pij = Ep[dWit.dWjt]/dt 
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T. By assuming that all the names share the same recovery rate R and with 
our previous notations, it means 

^pFtD{6) = (1 - R). - flil - 5, 

Invoking theorem |2] between t and t + dt and under the one factor copula 
model, the infinitesimal drift variation of a First- To-Default value is 

d^it,FtD ■■= _ 2\ X] \'^l^'it+f^%)PiPj-'^l^ttPitPij\ j Jl qk\x-(t>{di)(t>{dj)4>{x)dx, 

^ i,j,i<j k^i,i 

(14) 

if no default occurs in this time interval. Above, we recover the main formula 
of [5]. Assume that all the correlation levels are the same, in the real world 
and for pricing purpose: all the spread correlations pij are equal and are 
denoted by pg, and pi := p for all i. It is important to note that, whatever 
the values of correlations and spreads: 

• when p = 0, dpLt,FtD < 0. 

• when p = PS, dpt,FtD > 0. 

Thus, the break-even correlation of a First-To-Default lies somewhere be- 
tween zero and the corresponding (uniform) spread correlation. 

Actually, it is possible to find a general formula to deal with First p-th to 
default, p = 1, . . . , n. The default leg of such securities is the amount of loss 
in the portfolio, but capped at the p-th default. In other words, the payoff 
of a First p-th to default is 

n 

ipp = min(^ (5i,p).(l - R). 
i=l 

Such securities are important. Indeed, they may be seen as basic building 
blocks for most standard basket credit derivatives, including CDO tranches 
(see Brasch 2006). Now, we calculate the break-even correlations of First 
p-th to default explicitly. 

Theorem. 10 Under the assumptions (Cop) (I) (Q) (Fl) (A), the value of 
a First p-th to default will change between t and t + dt by the amount 

i<j 
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where 



and 



A;^j.7 a/1 -of, /I 



See the proof in the appendix |Ej Check that we recover equation ( 14 ) when 
p = 1. The previous reasoning we have led on the First- To-Default is avail- 
able for a general First p-th to default: in the case of identical spread corre- 
lations and identical pricing correlations, the break-even correlation level 
lies between zero and p'g. This result is true for arbitrary values of spread 
volatilities and every initial credit spread curves. 

At first glance, the theorem [TO] allows to write the Break-even correla- 
tion as a linear combination of spread correlations, with positive coefficients. 
Indeed, since all the quantities A*j above are positive, we can rewrite 

where all the Wijt belong to [0,1]. Nonetheless, the previous coefficients 
Wijt depend on pbe itself, through the so-called di terms. Moreover, the 
Wijt depend on the current values of the survival probabilities Qit and Qjt, 
i.e. they are time-dependent. Thus, the calculation of the instantaneous 
break-even correlation of a First p-th to default implies solving an implicit 
equation. 

If all the Qit are equal, i = 1, . . . ,n, then the previous terms A*j are 
equal. Thus, in this case, the (instantaneous) break-even correlation is given 
by 

2 

Pbe = ^ — j_ o2 \ ■ X] l^itl^jtPij- (15) 
2^i<j\Pit^ Pjt) i^j 

Note that this equation is true formally when there are two names only in the 
basket. In the latter equation, the break-even correlation depends explicitly 
on empirical correlations pij and some "volatility- type" coefficients Pi, that 
are equal to the previous quantities aiS^t (see section |3]). For instance, if all 
the spread volatilities and survival probabilities are the same, then 

2 _ 2 v-^ 

Pbe — ~( TT ' / , Pij ■ 
nin — 1 ^-^ 
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6 A bridge towards structural models 



All the previous analysis has been led inside the family of the copula based 
models. Traditionally, these models are defined through a direct specification 
of the joint laws of default times, but here, we have just needed the joint 
law of default events at the time horizon T. Now, it is fruitful to link 
this framework to structural models, following the seminal paper of Merton 
(1974). Classically, they are defined in terms of asset value processes (An) 
that follow some particular dynamics in continuous time. The default time 
of a name i is defined as a functional of z's asset value trajectory. 

The simplest case is given by the usual one-period Merton model, where 
in > T) := {AiT < biT), (16) 

for deterministic boundaries biT that are calibrated to be consistent with 
current observed default probabilities. Since the knowledge of asset values 
is equivalent of the knowledge of default likelihoods (or equivalently credit 
spreads), our asset values should depend on the same drivers as the survival 
probabilities Qu, i.e. on the previous Brownian motions (Wu), i = 1, . . . ,n. 
Typically, it is assumed the (risk-neutral) asset values follow some Ito pro- 
cesses 

dAit = Ti{t, Ait)dWit + Ui{t, Ait)dt, 

where (Wit) is a brownian motions under Q. For some classical families of 
diffusion processes the survival events can be rewritten 

iA,T < biT) =( [ i^^{u)dW,u < bir] , (17) 



for some deterministic function ipi and another default boundary biT- Note 
that the thresholds biT and biT, i = 1, ■ ■ ■ ,n, do not depend on t. In this 



section, we assume we can write (17). We deduce the survival probabilities 



that are implied by this structural approach: 

QMT) = Q(A.r < = 1. I ' I . (18) 



Now, in the previous Merton framework, assume that the correlation 
structure among the underlying names is one-factor: 



W^t = OiWt + \l - efW*^, i = l,...,n 



* geometric Brownian motion or Ornstein-Uhlenbeck process, e.g. 
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with independent Brownian motions (Wl-f.) under Q, and a common Q- 
Brownian motion (Wt) that generates some dependency between the un- 
derlying default times. Thus, the default events (rj < T) are independent 
given the r.v. Zi := Jg iJiiu) dWu/iJ^ i = l,...,n. If the r.v. Zi 

are not the same for different names i, the pricing procedure involves mul- 
tiple integrals, and can become tedious in practice. Hopefully, the usual 
Merton model is defined hy tpi = 1, and in this case, we recover the usual 
one-factor Gaussian copula specification. Otherwise, a Merton model with 
general functions ipi will induce m-factor models, m > 1 and m < n. 

Therefore, by defining default events through the relations 

{Ti>T)=( [ ^P^{u)dW^u < bir] , (19) 



\J0 / 

we get back the Gaussian copula model, possibly in a multi-factor version. 

Obviously, the available information at time t < T is the current sur- 
vival probabilities {Qit{T), . . . ,Qnt(T)), or equivalently the current asset 
values, or even {f^ ipilu) dWiu, ■ ■ ■ , Jq 'ipniu) dWnu)- Thus, we can write 
Tt = o (/g^ V'i(^) dWiu, i = 1, . . . ,n;v < ty Moreover, the common factor 
at time t (the "famous" X) is given by the trajectory of the common brow- 
nian motion (Wt) between t and T. To be specific, X can be reduced to the 
knowledge of the random variables ipi{u) dWu, i = 1, . . . ,n. 

Thus, let us define the filtration Q = {Gt)te[o,T]j Qt '■= c(/J" ipi{u) dWu, i = 
1, . . . ,n;v < t), and the extended filtration Ti = {Tt V Qt)t&[o,T]- 

For instance, in the standard gaussian copula model, the common factor 
X at time t is univariate and is simply Wt — Wt, the increment of the 
common brownian motion between t and T. In this case, ipi = 1 for every i. 
We prove in the appendix [F] 

Theorem. 11 In the IFGCM, the conditional probabilities {pj^^=w^-\Y^)t£[o^T]i 
as given in equation are %-martingale under Q. 

Now, we are interested in the Qjf-dynamics that are implied by the pre- 
vious Merton model. We would like to compare them with the dynamics we 
have been obtained in theorems [1] and [8j through fully different arguments. 



From (18), we deduce the dynamics of these survival probabilities under Q: 



dQit = -^it(^{^-HQit))dWu, Pit-= (20) 
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Note that is not random. Therefore, when to < t, we get easily 

to ^ to J t 

that tends to +00 when t tends to T. Actually, there is a one-to-one mapping 
between the deterministic functions f3i and ipi because 

Mt) ■= CI3it exp 1^ ^ du^ , (21) 

for some arbitrary positive constant C. 

Therefore, if the latter functions are of the type (3^ = d'iCt, or equivalently 



(22) 



then there is identity between the spread dynamics induced by the Merton 
approach (see (19) and (20)) and those obtained in theorem [s] by reverse 
engineering the GCM. 

For example, if the underlying asset values follow Brownian motions, 
then Ait ■= Wu and, with the previous notations, tpi{u) = 1 and biT = 
biT = VT^~^ {Qio{T)) . This is the most usual specification of a one-period 
Merton model. Then, in this case, I3it = Ct = l/\/T - t for all the names. 
With this very simple specification, the spread dynamics do not involve any 
firm-specific "volatility-type" coefficient a^: CTj = 1 for all i. 

Actually, the latter couple of functions {f3it,il^it) is not the single one. For 
instance, we could propose 

Wu,i^u) := (^^i^, - t/Tr^"^/'-'/'^ . (23) 

Note that this choice induces almost the same dynamics as in the previous 
case. But now, there is more heterogeneity in terms of spread dynamics, 
through different coefficients CTj. Therefore, to get heterogeneity (or "volatil- 
ity") effects, we can define the underlying survival events by 

in > T) := (^I\t - ut^^-^y^ dW,^ < 6.t) . 

Nonetheless, to price a basket in this model, it is necessary to integrate the 
conditional default probabilities with respect to several factors. 
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In the specification (23), it is nice to recover the results of section Wj 



With the notations of the latter section, the correlation between survival 
probabilities (or "spread" moves) pij is the correlation between the previous 
Brownian motions {Wit) Moreover, the correlations between the pricing 
factors were denoted by p'^pj when i ^ j r Since the asset values at maturity 
are given by Ait = Jq il^it dWu, i = 1, . . . , n, by definition of the vectors pi, 
we should satisfy p'^pj = Corr{AiT, Ajt)- But simple calculations provide 



Corr{AiT,AjT) 



Pij Jo i^iiJj 



Pij 



2ai(Tj 



(24) 



Thus, we recover the key relation of theorem |8j 

2a-,;a-,- 



PiPj 



Pij _ 2 I - 2 ' 



(25) 



Note that, through the equation (24), we have proved that, for an arbitrary 



correlation matrix [pij]-, every symmetrical matrix of the type 



Pij 



2(Tj(Tj 



is a correlation matrix. This was far from obvious. 
Another model would be to set 



(Ji 



(1 -t/T)"' {l-t/TY 



exp 



2(2a - 1 



-[{l-t/T) 



l-2a 



1] 



for some constants a < 0.5, ai > and ipi > 0. Under the latter Merton- 
type specifications, the asset value processes are different from usual Brow- 
nian motions, but follow (no drift, for instance) dAn = ipl^^dWu. Even in 
this case, we are still able to exhibit spread dynamics that allow perfect 
replication of CDO tranches. Through theorem [8j these are 

dQit = (3itUi^~\Qit))dW^t + pudt. (26) 
Simple calculations show that, again with the latter specification, we recover 



the usual relation (25) between spread correlations, pricing correlations and 
the "volatilities" a,-. 



see equation ( 20 



recall equation (111. The asset value processes are no more one-factor, but now p 



factors, p > 1, due to different functions Tpi 
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We have restricted our analysis to default events that are induced by 



structural models but through the relation (16), or, equivalently, through 



the relation (17). Other choices would be possible, for instance by defining 
a default event as a first hitting time of a boundary, i.e. (tj > T) = (A^ < 
biT, Vt G [0, T]). An open question would be to study such alternative struc- 
tural models, and to exhibit the spread dynamics that would be consistent 
with a "perfect" replication. Since multivariate extensions of First-hitting 
time models are very complicated in analytic terms (Zhou 2001, Hull et al. 
2006), this avenue is left for further research. 

It is interesting to note that the previous analysis has some striking 
applications in equity derivatives. A European digital put on z's stock has 
a payoff {SiT < K), for a given strike K and a maturity T. If we consider 
this stock follows a lognormal diffusion with constant volatilities cjj and a 
constant short rate r, then the t-price of this security is then 



:= QiSiT < K\Ft) = $ 



1 



ai^jT - t 



2 



ln(-|)-(r-^).(r-t) 



The dynamics of this price is the same as the classical one-period Merton 
model (when the asset values follow Brownian motions): 

dPit = -T^J i^'HPit)) dWu. 
\J 1 — t 

Moreover, the payoff of a digital worst-of European equity options is 
Vt4wo = '^{SiT < Ki, for some i G {1, . . . , n}). 

It is similar to the payoff of a First-to-default, after having defined the default 
events by {AiT < biT), i = 1, ■ ■ ■ ,n. Thus, the pricing formula of a First-to- 
default and the replication results we got in section [5] can be applied to the 
case of worst-of put options. For instance, at time t, the price of the worst-of 
put option on n names is 

Vt,,^o = 1-1 - ^) f^^7?4^) '^W (27) 

by assuming the correlations between stock processes are given by a p-factor 
model. Table [T] summarizes the correspondence between the credit deriva- 
tive and the equity derivative frameworks. If the true (historical) processes 
of stocks are lognormal diffusions, then historical correlations between these 
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Credit world 


Equity world 


Underlying 
Vanilla derivative 
Structured derivative 


asset value 
default probability 
First-to-default 


stock 
digital put 
worst-of put 



Table 1: One-to-one mapping between the European credit derivatives and 
the European equity derivatives. 



stock returns l^pij should be invoked in the pricing formula (27). Nonethe- 
less, it is not the case and the pricing formula (27) is still used, then these 
correlations should be different. For instance, if the stock returns follow 
diffusions 



d(ln Sj, 



(28) 



for some deterministic functions tl^i{t) = aiT{l — t/T)!'^^"^!/^ as above, then 
the pricing factors pi that will be put in formula (27) are now given by the 
relations 

, _ 2(Tj(Tj 



Nonetheless, a stock dynamics as in (28), even if not impossible, is unlikely 
because it depends on the option maturity T. Thus, a bootstrap-type pro- 
cedure would be necessary to build the whole dynamics of stocks. 



7 Alternative models 

In section |3] we have proved that the Gaussian copula model can be consid- 
ered (under some conditions) as a replication model if the individual survival 
probabilities follow some diffusion- type processes (theorem [2|. A natural 
question would be to find a similar result for other credit risk models, i.e. for 
alternative specifications of the quantities Surprisingly, it is absolutely 
not easy. Particularly, no Archimedean copula model can be seen as a repli- 
cation model, in our framework, as it is specified in the next theorem. Recall 
that a usual Archimedean (static) model is defined by a one-dimensional non 
negative random variable X, whose Laplace transform is denoted by ^, and 
by the conditional default probabilities 

^or equivalently between the put prices Pu 
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(29) 



for all > 0. 



Theorem. 12 For any Archimedean copula model (29), assume that 



• E[X^] < +00 



• the volatility functions an do not depend on i except through the tra- 
jectories Qiu, u <t. 

• The "real world" correlation levels pij at time t do not depend on the 
survival probabilities Qku, u <t, G {1, . . . , n}. 

Then there are no Qa- dynamics of type (Q) such that the price process {Vt) 
of any payoff ip is a T-martingale under Q. 

Actually, if we allow the spread correlations pij to be time-dependent, 
the picture is different. Indeed, some of the previous results can be extended 
to random correlation levels in a straightforward way. To be specific, we 
assume now that the i-instantaneous correlation level pij between dWu and 
dWjt can depend on the past realized values of the survival probabilities Qku^ 
A: G {1, . . . , n}, M < t. The main theorem [T] is still true because it is related 
to hedging strategies between two immediate successive times t and t + dt 
and conditional on past information 

Under this extended framework, we are able to exhibit an alternative to 
the Gaussian copula model, i.e. a "static" model that has an interpretation in 
terms of replication and the underlying single-name dynamics. The Clayton 
copula model provides such an alternative, even if it is an Archimedean 
model. 

We recall that, in the Clayton copula model, there exists a single common 
factor X that follows a Gamma law with the density fx{x) = x^^^~'^ exp(— x)/r(l/0) 
on for some nonnegative parameter 9. At time t, the default probability 
of i before T is given by 



In this model, the key dependence parameter is now 0, that is independent of 
i. Thus, the heterogeneity between names is driven here by the heterogeneity 
in terms of survival probabilities only. Now, we change slightly the dynamics 
of the Qit: 




(30) 
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Assumption (Q'): Under the risk neutral measure Q, the survival proba- 
bility process of name i is given by 



dQit{T) = audVu, t<T, (31) 
when i is not defaulted. Here {cru) are general J^-adapted process, and 

1/2 



dV^t := PiiQit) dZt + ^Jl- PHQit) dZ*t, pi{Q) := - (1 - Q) 

where Z and Z* , f = 1 , . . . , n denote n + 1 independent J^-Brownian motions 
under Q. 

The processes {Vit) are no more Brownian motions, in general. The 
following result is proved in the appendix. 

Theorem. 13 Under the Clayton copula model, assume (I) (Q') (Fl) (F2) 
(A). If 

o-jt = cro(l - Qit)- ( 1 - (1 - Qit) ) , i = l,...,n, 

then the pricing formulas given by the Clayton model are -martingales un- 
der Q. 

Thus, we have exhibited new spread dynamics that are consistent with 
the replication of CDO tranches, when they are priced with the Clayton 
copula. Moreover, we have linked the (unique) model parameter with the 
true spread process, through the instantaneous correlations between future 



spread moves. Note that, in theorem 13 when the default probabilities of 



name i goes to one, its spread correlations tend to zero. This is in line with 
intuition: all other things being equal, this name becomes more idiosyncratic. 

Moreover, under the pricing point of view, the previous Clayton model 
can be interpreted as a structural model one, like the Gaussian copula can 
be linked to the one-period Merton model. In the latter IFGCM case, the 
underlying asset values are driven by Brownian motions in the IFGCM. 
When dealing with the Clayton copula, we can consider asset values that are 
driven by a common Gamma process Y . Recall that, by definition, 

. y(o) = 0, 

• Y has independent increments, 

• Y{t) — Y{s) follows a Gamma distribution r(7.(t — s), A), t > s. 
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Recall that a random variable that follows a law r(7, A) has got the density 
X I— 7- A''' exp(— A2;)x'^~^/r(7) on Its Laplace transform is 

E[eM-sX)] = (l + ^ 

In |32) (rediscovered recently in |23|). a multivariate credit model can be 
defined by assuming the integrated hazard function of every name (say i) is 
given by 

Ai{t) = [ Ci{u)Y{du), 



for a given deterministic positive function a. In other words, at time 0, the 
default times are defined as 



Ti = mf{t >0\[ Ci{u) Y{du) > Ei}, 
Jo 



where Ei, . . . , En denote independent exponentially distributed random vari- 
ables. Equivalently, we can think in terms of asset values Xj, that would be 
defined by 



Y{t) 

Here, Ui, . . . ,Un denote independent random variables that are uniform on 
[0, 1]. In a one-period version of the associated structural model, the default 
event at time horizon T is defined by 

in <T) = {Xi{T) < h{T)), 

for some boundary bi(T). With this specification, the common factor is the 
random variable Y(T), and we can check that bi{T) = 1 — Qio{T). Thus, we 
deduce 

Pin < T\Y{T) = x,To) = exp (-xA[(l - Qio{T))-^/"'^ - 1]) . 

And we retrieve the equation (30) when A = 1 and 9jT = 1. 

In another perspective, we can try to find alternative replication models 
where every individual conditional probability is of the type 

Pi\x ■= V {^iX - HQit)) , (32) 

for some functions r], h and some constants Oj. The density of the common 
factor X is unknown, is denoted by fx as usual and its support is M. Is is 
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not a lack of generality: we could replace x by a function of a; in (32), and 
changing X's density conveniently. 

Note that the Gaussian copula model belongs to the latter family of 
specifications. It is the case of the Gamma-pool model too ([19], |22j). but 
not of Archimedean models. Thus, it is tempting to search some models 



of the type (32) that can be seen as replication models, i.e. that satisfy 
equation (j4]). Actually, it can be proved that the Gaussian copula only 
satisfies this requirement. 



Theorem. 14 Consider a model of the type (32), with unknown functions 
7], h and fx- Assume that the real world correlations pij are constant and 
that the volatilities an depend on i only through the survival probabilities 
Qi^,. This model satisfies equation ^ if and only if it is the Gaussian copula 
model, i.e. iff it is defined through equation 

See the proof in the appendix. 

Thus, the difficulty in finding alternative models can be seen as an ar- 
gument in favor of the Gaussian copula model. In other words, we can say 
that this specification has some nice analytical properties, that cannot be 
found elsewhere easily. The latter point should be stressed, when a lot of 
people criticize the Gaussian copula model (see Lipton and Rennie 2008, for 
instance). To summarize, the one-factor Gaussian copula model seems to be 
the most achieved one at least amongst of static models. Indeed, for this 
model only, we are able to 

• exhibit spread dynamics that are consistent with the model specifica- 
tion, 

• make an explicit link between the pricing parameters and the param- 
eters in the real- world dynamics, 

• price by replication when no jump-to-default occurs. 



To understand why the Gaussian copula model is so natural can be un- 
derstood by going back to structural models. If the underlying asset values 
Ait are driven by some Brownian motions, it is tempting to define the sur- 
vival of i at time T as the event "AiT is larger than some fixed constant". 
For most "natural" specifications, this is equivalent to "Wjt is larger than 
some fixed constant". In the latter case, we retrieve the Gaussian copula 
framework. To get alternatives, one way would be to leave the usual Brown- 
ian world, especially concerning the process that is followed by the survival 
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probabilities Qa. For instance, by assuming that the asset values (and thus 
the survival probabilities) are Levy-driven processes, but we have checked in 
theorem 14 that the Gamma-pool model does not satisfy equation ([s]), un- 



fortunately. Since these Levy-driven models allow spread jumps, the market 
efficiency can no more be obtained without considering alternative hedging 
strategies and hedging instruments. Another less natural way would be to 
define asset values still through usual diffusion processes, but by defining 
the default events as {AiT € li) for some unusual real subset Xj, or even as 
functions of the whole random trajectories {An), t G [0,T]. 

8 Empirical analysis of break-even correlations 

To illustrate these results, consider first a simple example: a First-to-default, 
second-to-default and third-to-default securities that are written on four 
names. We assume that the survival probabilities in the real world follow 
the "correct" theoretical dynamics ([8| without any drift. To simplify, we 
assume all the volatilities a are the same constant functions. Moreover, we 
assume that the dynamics associated to the first two names (say 1 and 2) are 
independent from the ones of the two other names (3 and 4). The "spread" 



between the two other names 



correlation between the first two names is = 30%, when it is p^^ = 70% 



It is relatively easy to check that the theoretical instantaneous break-even 
correlations for all the First p to default are given by 

2 _ ^12^*12 + ^34/^34 

Pbe 



A* -L. A* -L- A* -L- A* -L- A* 4- A* 
^12 ^ ^13 ^ ^14 ^ ^23 ^ ^24 ^ ^3 



with the notations of theorem 10 Obviously, the coefficients A*j^ depend on 
the order of the security, i.e. on p and on the current survival probabilities 
(or intensity levels) on every particular trajectory that is simulated. We 
draw some Qn, i = 1, . . . , 4, trajectories on 180 daily time intervals. Spread 
volatility coefficients at are assumed to be constant in time and the same per 
names. They take the value 50%. The ffist two names (resp. the last two 
names) share the same spot intensity level. Therefore, by averaging overs 100 
random trajectories, we are able to calculate mean break-even correlations 
on the whole time interval. 



If the spread trajectories are similar between all the names, then all the 
coefficients A*^^ are similar and the instantaneous break-even correlation is 
around (^^2 + pti)/^ ~ 16.7%. Note that this level is independent from 
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Spot Default intensities of 1 and 2 




0.1% 


1% 


5% 


30% 


0.1% 

■B ^ 1% 

1 s 5% 
-g CO 30% 


18 
45 
68 
69 


20 
17 
39 
69 


26 
17 

53 


30 

27 
18 
16 



Table 2: Break-even correlations of the First- To-Default in the four-name 
basket {a = 50%). 



Spot Default intensities of 1 and 2 




0.1% 


1% 


5% 


30% 


0.1% 

■B ^ 1% 

1 1 5% 
■g CO 30% 


18 
15 
11 
28 


12 
17 
15 
18 


15 
12 
16 
12 


32 
19 
11 
16 



Table 3: Break-even correlations of the First 2nd- To-Default in the four- 
name basket (ct = 50%). 



the strike we consider. Moreover, when the spreads of names 1 and 2 are a 
lot smaller than those of names 3 and 4 1^ then we can check easily on the 
theoretical expression ^ that 

• for the First- To-Default, « ^34 and ~ p§^ = 70%, and 

• for the Third-To-Default, « and ~ pfa = 30%. 

Indeed, these theoretical relations are observed in our simulation study: see 
tables |2] [3] and m 

Moreover, we have slightly modified the previous empirical framework: 
now, the basket has got ten names, and we have considered a particular 
trajectory of individual survival probabilities under ([8| (same time schedule). 
We have worked with constant volatilities at := (J for all t. We have assumed 
several scenarii in terms of individual spread volatilities, spread correlations 
Pij and spot intensity levels. Depending on the homogeneity among the ten 

*on the whole trajectory we consider 
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Spot Default intensities of 1 and 2 




0.1% 


1% 


5% 


30% 


0.1% 

■B ^ 1% 

1 s 5% 
-g CO 30% 


20 
17 
16 
29 


22 
18 
15 
25 


51 
25 
16 
19 


70 

69 
53 
18 



Table 4: Break-even correlations of the First 3rd-To-Default in the four-name 
basket {a = 50%). 



underlying names, we get different shapes in terms of break-even correlations 
(or equivalently, in terms of beta factors). 

The reference case is related to a fully homogeneous basket, with default 
intensities of Aq = 5%, constant volatilities of cr = 50% and constant spread 
correlations of pij = 50% for all i ^ j. In this case, the "break-even beta fac- 
tor" l^is closed to 50% and independent from the strike, as expected. Then, 
to move from this core scenario away, we have applied some deformations 
on these numbers: see tables |5] [6] and [7] Broadly speaking, these scenarii 
apply some rescaling method on all the parameters above, but differently 
from name to name. We observe different break-even correlation skews and 
smiles, depending on the heterogeneity of the underlying basket. For in- 
stance, in figure [T| different spread volatilities, spread correlations or initial 
default intensities per name (all other things being equal) generate down- 
ward sloping break-even correlations. This effect is more striking when high 
spread volatilities are associated with high spread correlations. Nonetheless, 
we can produce upward sloping break-even correlation skews by putting low 
beta factors on names that have got high intensities: figure [2] Note we get 
the opposite when the high beta factors are associated with high intensities. 
By combining these three dimensions, we can even generate some unexpected 
smiles/skews that are difficult to be forecasted: figure |3] 

We have calculated historical series of Break-even correlations, on two 
standard credit indices: ITraxx S8 Master, and CDX IG9 Main, both be- 
tween their launch in September 2007 and June 2009. In practice, we have 
calculated weekly P&Ls that are generated by continuously delta-hedging a 
base tranche [0, x%] (for ITraxx) or a senior tranche [x%, 100%] (for CDX). 
The hedging instruments are usual CDS on every name of the basket, and 

^i.e. the squared root of the break-even correlation 



29 



Name 


Core 


Up 


Down 


1 


0.50 


0.22 


1.13 


2 


0.50 


0.22 


1.13 


3 


0.50 


0.33 


0.75 


4 


0.50 


0.33 


0.75 


5 


0.50 


0.50 


0.50 


6 


0.50 


0.50 


0.50 


7 


0.50 


0.75 


0.33 


8 


0.50 


0.75 


0.33 


9 


0.50 


1.13 


0.22 


10 


0.50 


1.13 


0.22 



Table 5: Scenarii in terms of spread volatilities. 



Name 


Core 


Up 


Down 


1 


0.05 


0.02 


0.11 


2 


0.05 


0.02 


0.11 


3 


0.05 


0.03 


0.08 


4 


0.05 


0.03 


0.08 


5 


0.05 


0.05 


0.05 


6 


0.05 


0.05 


0.05 


7 


0.05 


0.08 


0.03 


8 


0.05 


0.08 


0.03 


9 


0.05 


0.11 


0.02 


10 


0.05 


0.11 


0.02 



Table 6: Scenarii in terms of default intensities. 
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Core 


Un 


T~)nwn 

I .y y J vv ±L 


1 


0.50 


0.22 


0.99 


2 


0.50 


0.22 


0.99 


3 


0.50 


0.33 


0.75 


4 


0.50 


0.33 


0.75 


5 


0.50 


0.50 


0.50 


6 


0.50 


0.50 


0.50 


7 


0.50 


0.75 


0.33 


8 


0.50 


0.75 


0.33 


9 


0.50 


0.99 


0.22 


10 


0.50 


0.99 


0.22 



Table 7: Scenarii in terms of beta factors. 



every relevant maturity. The CDO tranche calculation is led for a grid of 
constant correlation levels, from 0% to 95%. We calculate P&L increments 
with a rolling window of six weeks. For every window, the (empirical) break- 
even correlation is the level of flat correlation that allows a zero P&L vari- 
ation Moreover, to simplify, we assume that no interest are generated 
in the cash account, and that the running premia of all the instruments 
(CDS and tranches) are zero. Thus, they are managed fully with upfront 
paiements. 

For the sake of comparison, on our graphs, we have shown the break- 
even correlations and the associated base correlations, as deduced from the 



market: see figures |4| [5j |6j [T] for the CDX IG9 (main), and [8) |9j [lOj [ll| [12 
for the index ITraxx SB (main). It can be observed that the two series move 
consistently, even if the break-even series generate more volatile moves. But 
in average, the market base correlations are not so different from break-even 
ones, as expected, except for CDX senior tranches. It may be possible that 
investors in the latter tranches require an additional risk premium, or, in 
other words, overestimate systemic risk. 

In this practical experiment, several factors move us from the theoretical 
framework away: the CDO tranches are not European 

[3 heter ogeneity in 

terms of recovery rates for CDX, some default events have been observed 

^"Actually, to get nicer results visually, we have smoothed break-even correlations: at 
every date, the mapped break-even correlation is the weighted average of the current and 
three previous ones. The weights are exponential exp(— Afc), fc = 0, 1, 2, 3 and A = 0.3. 

^^thus, pricing formulas include immediate payments when default events occur 
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during this period of time for CDX , and have not been fully anticipated 



by the market. Moreover, our base correlations depend on the BNP-Paribas 



current analytics that have evolved in this period of time And, at the 
end of 2007, it has sometimes been difficult to find implied levels of base 
correlations from market quotes. Nonetheless, especially for ITraxx, it ap- 
pears that the concept of Break-even correlation can be a very useful relative 
value tool. Indeed, theoretically, its moves should be in line with those of the 
market base correlations. It is what we observe empirically, broadly speak- 
ing. Thus, significant and persistent gaps between both types of correlations 
should help to forecast future moves of base correlations, because both series 
should converge towards each other finally. 



9 Conclusion 

We have shown under which conditions usual "static" credit portfolio models 
can be seen as replication models. Explicit expressions of dynamic hedging 
errors have been provided for a large class of portfolio derivatives. In partic- 
ular, we have investigated the one factor Gaussian copula model, the market 
standard, and its multivariate extension. We have exhibited the unique fam- 
ily of dynamics (in terms of survival probabilities), that are consistent with a 
pricing by replication. We have shown that the matrix of spread correlations 
Pij and some spread "volatility-type" coefficients ai should be combined in a 
particular way to get break-even correlations. Broadly speaking, the matrix 
[2(JiajPij/{af + a'j)] should be definite positive. When it is the case, this 
matrix is our break-even correlation matrix and should be used for pricing 
and hedging purpose. We have shown that the one-factor copula is more 
tractable than most other factor models and why it is so difficult to find an 
alternative to it by proving that most other models will not lend them- 
selves to replication. 

There remains a lot of avenues for further research. In particular, we 
should try to adapt our framework to deal with 

(i) richer dynamics in terms of survival probabilities, for instance by as- 
suming underlying Levy processes, 

(ii) random recoveries (Amraoui and Hitier 2008, Krekel 2008), 
^^Fannie Mae, Freddy Mac, Washington Mutual 

^^particularly, a switch towards a stochastic recovery model in the last quarter of 2008 
^*at least inside the family of static models 
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(iii) some extensions of the one-factor Gaussian copula framework, as ran- 
dom factor loadings (Andersen and Sidenius 2004) or local correlation 
(Turc et al. 2005), 

(iv) perfect hedging when sudden jumps-to-default may occur before ma- 
turity. 

(v) joint specification of all the Qjt(T)-dynamics for different maturities 
T, i.e. an arbitrage-free specification of term structures. 

All these points induce significant technical difficulties. Concerning the point 
(i), additional jumps components would challenge the desired perfect repli- 
cation property. Moreover, random recoveries cannot be tackled easily. At 
ffist glance, it could be done by working with individual loss dynamics in- 
stead of survival probabilities Qit{T). In point (iii), the idea is to build a 
dependence between the pricing correlation level and the common factor, 
i.e. to price with some functions Pi{x), i = l,...,n. Point (iv) has been 
discussed in section [2j Since it should involve alternative sensitivities and 
alternative instruments, this task has been left for future research. The last 
point (v) is probably the most crucial one. All our analysis has been led 
with a given maturity T. With two maturities T < T' and the same frame- 
work, we can check that it is not possible to satisfy the arbitrage relation 
Qit{T) > Qit{T'), for all times t < T. Indeed, when t tends to T, Qit(T) will 
tend to zero with a certain probability, when Qit{T') will not exhibit such a 
feature (as a consequence of the spread dynamics we have found). To find 
other model specifications, similar to ffist-hitting time models for instance, 
is the main challenge to tackle the multi-period formulation of the problem. 

The authors thank their colleagues and former colleagues of BNP-Paribas 
and JP-Morgan for their help and fruitful discussions. Bruno Bouchard, 
Stephane Crepey, Peter Jackel, Monique Jeanbleanc, Jean-Paul Laurent, 
Dong Li and Marek Rutkowski have provided us highly valuable comments 

too. 
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A Proof of theorem [l] 

By the conditional independence property, the price of the structured prod- 



uct is 




with our notations. We have 



d'Vt 




and, for every couple {i,j), i ^ j 



d^Vt 



dQi dQj 



X 



fx{x) dx. 



dQidQ 



(34) 
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Then, by applying Ito's formula, we get 

dVt = {...)■ dWt + 



2 ^ dQidQj 



r 2 f)2 

5{T) ^^^hj * 

:= (...)• dWt + (iTTi + dvrs. (35) 
By an integration by parts and with our notations, we get for every i 

r n -l-oo 

We deduce 

and the formula follows. □ 





+00 






— oo 





B Proof of equation ([6j) 

To lighten the notations, set Sj := ^~^{Qjt)- Simple calculations provide 

dPj\x _ (pjdj) 



and 

d'^Pj\x _ (pidj) 



3' 

Thus, 



Xjix) = „ — 2v?72 / ~ P3Sj\.ct>{dj)(t>{u) du. 



(,.)2(l_p2)3/2 
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But a key point in favor of the Gaussian copula specification is that 



^{dj)cPix)=<P\^7==\<f>i~^j)- (36) 



Therefore, we deduce 



_ pj(j){dj)(f){x) _ Pj(i){x) dpj\x 
Moreover, we get easily 

-Pj(l){^ HQjt))- 



and the result follows from theorem [H □ 



C Proof of theorem [H 

Clearly, conditions (i), {ii) and [Hi) are sufficient so that iVt) is a J-- 
martingale under Q, invoking theorem [2] 

Now, let us prove the necessity of such conditions. We want to cancel the 
drift of dVt for any random trajectory of all {Qit)tel, I C [0, T), z = 1, . . . , n. 
To fix the ideas and without a lack of generality, we consider the particular 
couple of indices (1,2). In particular, the drift of dVt should be zero for 
trajectories where the Qa take very large or very small values, i 7^ 1,2. To 
be specific, we will consider trajectories such that sup^, (j){di)(j){x) < e, for a 
given e > 0, and when i ^ 1,2. Thus, since we are under the "no default" 
assumption, under our assumptions, the drift of dVt is 

I [2l3ul32tPi2 - PiP2l3ft - P1P2PI] J2 ^i^TmSi - l){252 - 1) 

5{T) 
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We have assumed that there exist some indicators (5^, k ^ 1,2 such that 

^'(i, 1, <ri(,^2)) - v(i, 0, ^"'1(1,2)) - v(o, 1, 51(1^2)) + ^'(o, 0, r (1^2)) + 0. 

Now, consider a compact subset A such that J l(x ^)(/)(x) do; < e. We 
can particularize our trajectories Qa even more, so that 

sup J] < 

when (5fc / (5^ for at least one index /c / 1, 2, and 

This is possible, because, for every index A; 7^ 1, 2, we can find Q\^i such that 
p^l^ (resp. is closed to one when b\ = \ [b^. = 0), and uniformly w.r.t. 

a; in a compact subset. 

Therefore, the sum defining the hedging error is "reduced" even more: for 
this particular choice of Qkti ^ 7^ hJi the drift of dVt is 

I (V(l, 1, 5"** (1,2)) - i^ih 0, 5l(i,2)) - ^(0, 1, 5"** (1,2)) + m 0, ^1(1,2))) 

ro« « «2 fl2 1 /■ TT ^fc (/)(di)(?:)((i2)</'(x) , N 

• [2/3ii/32tPi2 - PiP2^u - Pip2p2t\ / 11 Pkix'iklx r, 2 rr^ ^ + '^^^^ 

fc^l,2 ~ PlV^ ~ P2 

= Cst [2/3it/32tpi2 - pmPlt - Pip2^lt] / (t){di)(l){d2)(t){x) dx + 0(e). 
Since e is arbitrarily small, this drift is zero if and only if 

2PuhtPl2- PlP2[Plt + Plt\=^. 

for every Qu and Q2t- 

Actually, it is possible to lead the same reasoning for all couples (i,j), 
i ^ under the assumption (P). Thus, we recover equation Q, that is now 
a necessary and sufficient condition: 

2PitP,tPij = [Pl + Pl]piPj, (37) 

for every x, Qjj, Qjt and every couple (i,j), i 7^ j. The latter condition 
can be true if and only if the previous (in are independent of Qu- They 
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can depend on time t, but iff Pit/f3jt is independent of t for all (i, j), i 7^ j. 
Therefore, we get equation (|8]), by rewriting f3it as aj^f, for some positive 
function ^ and some constant > 0. Our assumption on S"^ provides the 
existence of the pricing parameters pi, i = 1, . . . , n. 

To finish the proof, it remains to prove that lim^-i.^ /q du = +00. By 
setting Zf = ^~^{Qit) and invoking Ito's lemma, we see easily that the 
solution of the PDE ([8| is given by 

/ f* f 2 du r* r* ^2 \ 

Qu = $ ( exp(^^^o^)cI>-i(Q,o) + eM ' P^CudW.A , 

where (Wiu) is a Brownian motion under Q. Note that Qit{T) tends to zero 
or one when t tends to T. Assume that limt^T foCud-u = i < +00. Then, 
exp(af ds/2)S^udWiu is a non degenerate Gaussian r.v. and it will 

be impossible to get liuit^T Qit G {0, 1} almost surely. So the result. □ 



D Proof of theorem [6l 

We come back to the first steps of the proof of theorem [T] in section |A] With 
these previous notations, we have now 

and 

(5(r) '^^^.i * 

(39) 

Under the p-factor Gaussian copula framework, these quantities can be 
rewritten: 

^(y) i=l i^i ^ l'^*! ^ fc=l 

(40) 



where, as previously, we have set Si := $ ^{Qit)- Moreover, 
d7T2 _ 1 

~dt ~ 2 



E E / n 4Urx '=•(2'^^ - 1)(25. - 1) 

(5(T) 

^{di)<^{dj) nLi '^(^fc) ^^fc /..N 

• "^"^■^^^>(.-.)0(.-,) ■ v^r^wyr^^' ^ ^ 
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The key point of the proof is to exhibit the right variables to make in- 
tegration by parts. Let us introduce some additional notations : for any 
vector X in and r S {1, . . . ,p}, let us denote by x_,. the p — 1 vector 
{xi, . . . ,Xr-i,Xr+i, ■ ■ ■ ,Xp). For every integer r < p, simple calculations 
provide the identity 



(42) 

where Si^r = Si — p[ _^x-r- Now, it can be checked that we can rewrite 



p 

P'i^ ~ \Pi\'^Si = [pirXr{l — \Pt-r\ J — ^ir \ j 



where we have set 



2 - 2 / 

Cir •= Pir^i ~ PirPi—r^—r- 

We deduce dvri = X]r=i dt^i^r-, with 

d-Kl^r sr^ irT^rr^w-^ [tJ <5) I--?? /o'T 'Pidi)[PirXr{l - \Pi-r\'^) - Cir] x , 

^ = E^(W)E y n^MixM2<^.-i)f — ^(,-^)2(i _ 1^ ,2)3/2 — ii^(-^)dx,. 

5(71) 1=1 jy^i T \ / \ ir I / J.^;^ 



But, by integrating w.r.t. the variable Xr and by invoking equation (42), we 
observe that 

4>{di)[pirUr{l - \pi-r\'^) - Cir]4>{Xr) dUr 



/I I |2\ / ( / IPij-J'P r Pi,rSi r 



Pir {l-\pi\^)ct){d,)(P{Xr). 
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Thus, by an integration by parts w.r.t. x,-, we get 



5{T) *=i '^^''i 

/9i^(l - \pi\'^)(l){di) -A- 

• .2/1 1 12^3/2 11 '^(^fc)^^fe 

'/'(•Si) (1 - j'^/ 

= - E V'(w) E / n 4Urx ^-(^-^^ - - 1) 

S{T) 

PirPjr(t>{di)4>{dj) 

0(.-)2(l-|p,|2)l/2(l_|p^.|2)l/ 



And then, by summing up on r, we obtain 



5{T) *<i '^T^^.i 



By summing the equations (43) and (41) up, we get the result. □ 



E Proof of theorem 10 



Let i;p{6) = mm{J2'l^^6up).{l - R), and 6p{6) = ^E^=l^^ < P)-i^ " 
We can checked easily that, for every p > 1, 

ijjp = ipp-i + 1- R- 6p-i. 

The drift of the value process associated with a particular payoff ■0 will be 
denoted by d'K{ip). Obviously, vr is a linear functional, and d7r(V') = if the 
payoff is risk- free. We deduce 

dn{^p) = dTT{iPp^i) - dTT{9p.i). (44) 

Now, let us specify dTi{Op) for p < n. To lighten the notations, we set 

^ij - ^ [iPitPjtPij - PiPjPit - PiPjPjt] 



42 



Let us apply theorem [2] By distinguishing between the different values that 
are taken by the couples {5i,6j), we get 



5(T) '=1 k^iJ 



EE 



l(5^<5/<p-2)-21(^5i <p-l) + l(j;5, <p) 

l¥'id l¥^ij Ij^ij 



where denotes the vector of {6k), with k £ {1, . . . ,n} but k / To 

simplify, let us denote 



Oir. 



If p > 1, we have got 

dMOp) = Yl E [i(E'^'=^)-i(E'^'=^'-i)] / Ilpl\A\'J''^,,dxdt, 



We deduce, for all p > 2, 

dTr{ipp) = d-K{il)p-i) - ap-i + Qp_2- 

By iterating the latter equation, we get d'K{tpp) = d'K{tpi) — Up-i + oq- But 
7r(V'i) is the hedging error of the First To Default, that is equal to — (see 



equation (14)). Finally, we get dir^ipp) = —Op-i, so the result. □ 



F Proof of theorem [TT] 

In the one-period Merton model, with our notations and given Tit, 



pj,t\x = Q (r, < T\nt) = Qyj^ ^M) dWju < bjTint 

f-T ft 



Q[Jl-p'. / ^j{u) dW*^ < bjT - / i^j{u) dWju - Pj / Mu) dWM 







^ I - lo i^j ju) dWju - Pj Jt i^j (n) dWu . _^ 
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When ipj{u) = 1 for any j and u, we recover the standard Merton model and 
the pricing formula (|5]). In this case, Ht = (T{Wju, Wt — Wu, u < t). For any 
t,t', T >t > t' > 0, we get 



E 
E 



1 - pjVr^tz + {w. 



jt 



jt') 



1 - pjVr^tZ +Jl- pjiW*t - W*t,) < hr - Wjt' - PjiWr - Wf)) \n 



for some gaussian r.v. Z, independent of (Wjt) and (Wt). We get 



E\Pj,t\x\'Ht'] = $ 



bjT-Ww -Pj{WT-Wt' 



'l-pjVT^ 

Therefore, the process {Pjt\x=WT-Wt)t(^[o,T] is an 7^-martingale under Q. □ 



G Proof of theorem [12 



In an Archimedean framework, we would like to find some dynamics that 
allow a cancellation of the hedging error dvrj over all trajectories (Qu), 
i = 1, . . . ,n, t G [0,T] and for all payoff -0. It implies we have to satisfy 
equation Q by the same arguments as in the proof of theorem [C] Partic- 
ularly, we assume the case Qit = Qjt for all t and for some couple 
^ 7^ J E3 Since the volatilities are the same measurable functions of their 
trajectories (Qu), this implies that au = (Jjt on these particular trajectories. 
Thus, we have to satisfy the necessary relation 



dp, 



'i\x 



dQi 



fxix) = Xi{x) 



dp, 



'i\x 



dx 



(45) 



for some non zero constant a. By simple calculations, it is equivalent to 

l2 



ax'^Pi\xfx{x) = [^'{Ci)] ^iXi{x) 



^^Even if this event occurs with probability zero, the latter assumption can be justified. 
It is sufficient to consider trajectories s.t. \\Qit ~Qjt\\ < £ for some arbitrarily small e > 0. 
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for all X in the support of X and every Qi G [0, 1]. By deriving w.r.t x and 
by making the change of variable y = we have to satisfy 

a[2-yx + x^ix)]=y[x + ^{y)], (46) 

for all X in the support of fx and for all y G [0, 1]. Particularly, when y tends 
to zero, jr{y) tends to — £^[X^]/£'[X] that is finite by assumption. Thus, we 
have to satisfy 

f 

x^(x) + 2 = 0, 
Jx 

for all X > 0. This implies the support of the density function fx contains 
some interval [a, f3] C M+, where it is equal to 

fx{x) = 



for some positive constant a. But then, equation (46) implies a = ( — 1) and 
= 0, that cannot be satisfied clearly. □ 



H Proof of theorem 13 



Under our assumptions, we calculate 

-xOpiixCi- - Qit) 



d Pi\x 



{xeypi\x{i - Qit)-''-' - xe{e + i>i|,(i - q 
%^=P.|.[i-(i-Q.*)-% 

By an integration by part, we obtain the identity 

dPi\x x9fx{x) 



it 



Xi{x) 



dQi 1 - Qi 



Therefore, with our choices of volatilities an and the underlying processes, 
it can be checked easily we satisfy (Q, replacing pij by 

< dQit,dQjt >= critcrjtf3'i{Qit)/3j{Qjt). 

Finally, we do not have to take into account the so-called term r] of the 
theorem [T] because liuix^o Xi{x) = and limi._>+oo Xi{x) = 0. □ 
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I Proof of theorem 



14 



To lighten the notations, we set rji := rj {atx — h{Qit)), hi := h{Qit) and 
/ := fx{x). Under our assumptions, we get easily 



dpi\x 
dx 



dQi 



-n'h' ^ 



vA- 



We would like to satisfy the partial differential equation Q that is rewritten 
here 

(47) 



2{aih'i).{ajh'j)pijf = ajaf^ + aia]^ 



2X3 



The latter equation has to be satisfied for all the triplets {Qi,Qj,x). Actu- 
ally, it is a very strong constraint, because it can be proved that this implies: 
there exist some constants Cj and C{i,j), i, j = 1, . . . ,n, such that 



CTih'i 
2Xi 



Ci, 

C{i,j)pijf, 



(48) 
(49) 



for all indices (i, j), i 7^ j- The latter result is a consequence of the following 
lemma: 

Lemma. 15 Consider some real functions a, b, c, d such that 

a{x)h{y) = c{x) + d{y), 
for all couples {x,y). Then these four functions are constant. 

This lemma is proved at the end of the current proof. Note that d is not 
a function of x because ah'^ is a function of Qi only (not of x). Moreover, 
equation (49) is due to the fact we can consider the particular case ai = aj 
and some trajectories where Qn = Qjt 

From the latter equations, we get the existence of a constant p = p{i,j) 
such that 

Xi = f^fViih'i)^- 
By differentiating this equation w.r.t. x, we get 



M2 



+ ai- 



(50) 



^In this case, the two members of the rhs of equation (481 are equal. 
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Now, by integrating the latter equation w.r.t. x, we obtain 



r?^ = exp ( + ) (51) 

(1 -^ai)(/ij^ 



for some function i/j := v{Qi). But (50) can be multiplied by h'^. And after 
an integration w.r.t. Qj, we get 

.;^e.p(;M- t'"W+M.w/] y ,,,, 

V (1 - ^iai) J 

for some function S,- If /'// were a constant function, the density / would 
be exponentially distributed and the support of X cannot be M as a whole. 



Moreover, /i is a not a constant by assumption. Thus, equations (51 ) and (52 ) 
together imply that the crossed terms f'hi/f and xh'-(h[)~'^ are proportional. 
In other words, this means that /'// (resp. hih[) is proportional to x (resp. 
h'-/h'-). Therefore, by integration, we check that the usual Gaussian copula 
model is the single possible model that satisfies such requirements. □ 

To finalize this proof, it remains to state lemma [TSj This is done now. 



If there exists some real number x such that a{x) = 0, then it is clear 
that c and d are constant functions, and then a and b too. We can lead the 
same reasoning with b. Thus, we can assume that a and b cannot be zero. 
Now, assume that b is not a constant: there exist a couple (yi, 2/2) such that 
b{yi) 7^ b{y2) 7^ 0. Then, for every real number x, we have 

Kvi) ^ c(x) + d{yi) 
Ky2) c{x)+d{y2y 

that is constant function in x. This implies that c is a constant function 
of X, and that d = ab + 13 for some constants a and (3. We deduce that 
a{x) = (c + ab{y) + f3)/b{y), so a and b are constant functions. 

This is a contradiction with our initial assumption. Since this reasoning 
can be lead with a instead of b similarly, we have proved our lemma. □ 
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Figure 1: Break-even "beta factors" on a simulated trajectory 
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Figure 2: Break-even "beta factors" on a simulated trajectory 
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Breakeven Correlation Skews/ Smiles 




1 2 3 4 5 5 7 8 9 

Detachment point 



— Fully homogeneous — High vol, low I ntensity, high spread correl — High vol, low intensity, low spread correl 



Figure 3: Break-even "beta factors" on a simulated trajectory 
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Break-even and Base correlations: CDX IG9 3-100% 




25-Sep-07 25-Dec-07 25-Mar-08 25-Jun-08 25-Sep-08 25-Dec-08 25-Mar-09 



—Break-even correlation —Base correlation 



Figure 4: Break-even correlations and base correlations, CDX IG9, 3% — 
100%. 
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Break-even and Base correlations: CDX IG9 7-100% 




25-Sep-07 25-Dec-07 25-Mar-08 25-Jun-08 25-Sep-08 25-Dec-08 25-Mar-09 



—Break-even correlation — Base correlation 



Figure 5: Break-even correlations and base correlations, CDX IG9, 7% — 
100%. 



52 



Break-even and Base correlations: CDX IG9 10-100% 
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Figure 6: Break-even correlations and base correlations, CDX IG9, 10% — 
100%. 
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Break-even and Base correlations: CDX IG9 15-100% 
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Figure 7: Break-even correlations and base correlations, CDX IG9, 15% — 
100%. 
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Break-even and Base correlations: Itraxx S8 0-3% 
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Figure 8: Break-even correlations and base correlations, ITraxx Main S8, 
- 3%. 



55 



Break-even and Base correlations: Itraxx S8 0-6% 
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Figure 9: Break-even correlations and base correlations, ITraxx Main S8, 
- 6%. 
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Break-even and Base correlations: Itraxx S8 0-9% 
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Figure 10: Break-even correlations and base correlations, ITraxx Main S8, 
- 9%. 
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Break-even and Base correlations: Itraxx S8 0-12% 
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Figure 11: Break-even correlations and base correlations, ITraxx Main S8, 
- 12%. 
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Break-even and Base correlations: Itraxx S8 0-22% 
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Figure 12: Break-even correlations and base correlations, ITraxx Main S8, 
- 22%. 
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